Viscoelasticity and primitive path analysis of entangled polymer liquids: 

From f-actin to polyethylene 
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We combine computer simulations and scaling arguments to develop a unified view of polymer 
entanglement based on the primitive path analysis (PPA) of the microscopic topological state. Our 
results agree with experimentally measured plateau moduli for three different polymer classes over 
a wide range of reduced polymer densities: (i) semi-dilute theta solutions of synthetic polymers, (ii) 
the corresponding dense melts above the glass transition or crystallization temperature, and (iii) 
solutions of semi-flexible (bio)polymers such as f-actin or suspensions of rodlike viruses. Together 
these systems cover the entire range from loosely to tightly entangled polymers. In particular, we 
argue that the primitive path analysis renormalizes a loosely to a tightly entangled system and 
provide a new explanation of the successful Lin-Noolandi packing conjecture for polymer melts. 

PACS numbers: 61.25.H-, 83.10.Kn, 82.35.Pq, 87.10. Rt, 87.10.Tf 



I. INTRODUCTION 

The relation between the complex viscoelastic proper- 
ties of polymer liquids and their microscopic structure 
and dynamics is a key issue in materials science and bio- 
physics [ll, 0, H, 0, On a microscopic scale chains 
can sHde past each other, but their backbones cannot 
cross; the Brownian motion of these macromolecules is 
hence subject to transient topological constraints Q, an 
effect which is familiar from the manipulation of knotted 
strings. In slowing down the chain equilibration after 
a deformation, these constraints or entanglements dom- 
inate the viscoelastic behavior of high molecular weight 
polymeric Hquids. Entanglement effects are universal, i.e. 
one observes the same behavior for polymers with similar 
overall chain architecture (linear, ring, branched) inde- 
pendently of details of the molecular structure. Modern 
theories of polymer dynamics and rheology [J, [3| describe 
the universal aspects of the viscoelastic behavior based 
on the idea that molecular entanglements confine indi- 
vidual filaments to a one-dimensional, diffusive dynam- 
ics (reptation) in tube-like regions in space. Material 
specific parameters are determined through comparison 
to experiment. Here we are concerned with the question, 
if these parameters and related experimental observables 
can be inferred from the molecular structure of polymeric 
liquids. 

How strongly linear polymers entangle with each other 
depends on their stiffness and on the contour length den- 
sity of the polymer melt or solution . The microscopic 
structure is best discussed in terms of the Kuhn length, 
Ik, and the number density of Kuhn segments, pK- The 
Kuhn length is defined as the contour length L where 
thermal fiuctuations start to bend the chains and marks 
the crossover from rigid rod to random coil behavior. In 
"loosely" entangled systems with pkI\c < 1 the mean- 



free chain length between collisions is larger than the 
Kuhn length, leading to random coil behavior between 
entanglement points. In contrast, for pkIk ^ 1 ^1^" 
ments are "tightly" @| entangled and exhibit only small 
bending fiuctuations between entanglement points. As a 
consequence, the chains in tightly entangled f-actin solu- 
tions and in loosely entangled polyethylene melts behave 
differently on the tube scale. The former are essentially 
stiff and resist macroscopic shear due to an increase of 
their bending energy. The latter are fiexible and lose 
entropy when stretched locally. 

The differences between these two situations become 
apparent, when one considers the relation between the 
microscopic solution structure and the height of the 
characteristic rubber-like plateau in the shear relaxation 
modulus, G%. In Fig. [l] we show a Graesslev-Edwards 
plot 01 of experimentally measured 0, [13, [U, H [13, Q 
dimensionless plateau moduli G^l^/kBT as a function 
of the dimensionless number density of Kuhn segments 
PkIk- comparison, the unsealed plateau moduli are 
shown in the inset as a function of polymer concentra- 
tion. The data we have compiled represent the behav- 
ior of several prototypical classes of entangled polymers: 

(i) tightly entangled solutions of semi-fiexible biopoly- 
mers such as f-actin and suspension of fd-phages with 
PkIk > 10, « IQ-^m and G% « lO-^Pa, and 

(ii) loosely entangled melts of commercially important, 



synthetic polymers with pkIk < 10, 



10- 



and 



plateau moduli of the order of G^^ ~ lO^Pa. Typical 
synthetic polymers are sufiiciently fiexible to be in the 
isotropic phase in the melt state with a volume fraction 
^ = PK Ik d'^ — 0{1). In contrast, tightly entangled 
chains have to have a sufficiently small diameter d to ful- 
fill the Onsager criterion ^{iK/d) < 5 for the isotropic-to- 
nematic transition. We note that the crossover between 
the two regimes in Fig. [J is located close to the thresh- 
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FIG. 1: Dimensionless plateau moduli G% Ix/ksT as a func- 
tion of the dimensionless number density of Kuhn segments 
PkIk- Experimental data vs. scaling predictions. Symbols 
indicate experimental data for (i) ti ght ly entangled f-actin 
(blue +) [a] and fd-phage (blue x) [1CI | solutions, (ii) var- 
ious loosely entangled polydiene, polyolefine, and polyacry- 
late melts (green x) [ll|, u^, and loosely entangled theta- 
solutions of polystyrene (red -I-) [l^ and polybutadiene(red 
x) Lines represent theoretically deriv ed p ower laws 

G% oc [pkIk)" for tightly [3| and loosely |l6l | entangled 
systems. 



old for the isotropic-nematic transition in dense polymer 
melts [IBl- Furthermore, Fig. [l] contains data for (iii) 
loosely entangled semi-dilute solutions of synthetic poly- 
mers in so called theta-solvents. While these systems are 
supposed to preserve the chain conformational statistics 
from the undiluted melt, the observed reduced plateau 
moduli are larger than those for dense systems and ex- 
hibit a qualitatively different density dependence. 

For tightly entangled systems, the relation between 
the entanglement density, pe, and pkIk determined 
by Semenov jl^ from a geometrical argument: the area 
swept out via transverse fluctuations by a fllament be- 
tween two entanglement points is on average traversed 
by one other filament serving as an obstacle. In loosely 
entangled systems it is impossible to determine Pp by 
scaling arguments alone 0, M, M, [13, E E H MMM- 
A promising tool to solve this problem from first prin- 
ciples is the primitive path analysis (PPA) [H, [23, H, 
m, HE, [131 ■ Primitive paths were originally introduced 
in a thought experiment to determine the tubes confin- 
ing individual polymers in an entangled polymer melt 
or network. The idea is to identify the random walk- 
like tube axis with the shortest ("primitive") paths be- 
tween the end points of the original chains into which the 
chain contours can be contracted without crossing each 
other jH, [l^l- As we have shown [.1.3, 23], the numerical 
implementation of this idea allows to make quantitative 
predictions of melt viscoelastic properties on the basis of 



a topological analysis. Here we (i) argue that the PPA 
can be understood as a way of renormalizing a loosely 
to a tightly entangled system, (ii) derive a modified rela- 
tion between the primitive path mesh characterizing an 
entangled polymer liquid and its macroscopic properties 
(Sec. [ni), (iii) apply the PPA to a much wider range of 
loosely and tightly entangled model polymer structures 
(Sec. imp , (iv) validate our results by a comparison to the 
experimental data displayed in Fig.[ll and (v) discuss the 
relation between the characteristic length scales of the 
primitive path mesh and the packing length (Sec. lIVp . 
We conclude with a brief summary in Sec. |Vl 



II. THEORY 

The result of the PPA is a mesh of mutually entangled, 
piecewise straight primitive paths (see, e.g. Fig. 3 in 
Ref. [131 )■ The structure can be characterized by the 
contour length, Lpp, the mesh size, £,pp — 1/ \/ PchainLm, 
and the Kuhn length, Opp, of the primitive paths [l^.[23|. 
Furthermore, it is useful to introduce the average contour 
length, Ze, of the primitive paths between entanglement 
points as well as the corresponding chain contour length, 
Lg. On a scaling level, the latter is defined implicitly by 
the relation 

11 = {R^){L,,Ik) (1) 

where {R^){L,Ik) denotes the mean-square spatial dis- 
tance of two points separated by a length L along the 
contour. The chain and primitive path statistics agree 
asymptotically limL->oo(^^) = IrL = appLpp — 
limLpp^oo(^^)pp so that the Kuhn length 

Ik {le/Le) 

of the primitive paths increases by the inverse of the 
shrinking factor of the contour length. Similarly the 
primitive path mesh size is given by 

Ik {pKiiy^ielL^yn- ^' 

In the following, we assume (i) that all information 
necessary to calculate the plateau modulus of a sample 
can be deduced from the primitive path mesh character- 
izing its microscopic topological state and (ii) that the 
primitive path structure - viscoelastic property relation 
is system independent, i.e. after carrying out the PPA 
it is no longer necessary to distinguish between different 
polymer classes. Moreover, we note that by construction 
primitive path meshes resemble tightly entangled solu- 
tions of semi-fiexible chains and that the latter are in- 
variant under the PPA. 
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Regarding the primitive path analysis as a means to 
renormahze a loosely to a tightly entangled system, al- 
lows us to adapt two results from the theory of tightly 
entangled solutions of semi-flexible chains |15| to the 
present situation: 



'^N '■K 



C5 le X {ll/appf^ (4) 

CGcl'\pKliy"'{lelL,f"' (5) 



Eq. (U expresses the idea that the area swept out via 
transverse fluctuations by a primitive path between two 
entanglement points is on average traversed by one other 
primitive path serving as an obstacle. Eq. (0 states that 
the plateau modulus is proportional to fc^T times the 
density of entanglement points. 

For given chain statistics it is possible to determine le 
self-consistently from Eqs. |(T]|4l). For tightly entangled 
systems, our approach reduces by construction to the 
standard results for semi-dilute solutions of semi-flexible 
chains [H, (s^]: app/lK = 1, ^.pp/lK = 1/{pkIkV^^, 
yiK cx and Glll/ikBT) cx {pkI'kY^'- In 

the opposite limit, we flnd Upp/lx cx: £,pp/lK oc 1/{pkIk) 
and G%l\/{kBT) cx {pk^k)^ i-^-) we have presented a 
derivation of the Lin-Noolandi conjecture [l3,[23|- In the 
general case of solutions where the individual polymers 
exhibit worm- like chain statistics on all length scales, a 
reasonable approximation (indicated by solid lines in all 
flgures) is given by 
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Ik 



Le 

Ik 
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-1/2 



c^PkIk) 



2/5 



c^PkIk) 



(6) 
(7) 



By inserting Eqs. © and into Eqs. 0, P, and © 
we thus arrive at a prediction for the dependence of the 
properties of the primitive path mesh on the dimension- 
less Kuhn length density pkIk of the original polymer 
melt or solution. 

Choosing = 0.06 and cg = 0.6 we flnd excel- 
lent agreement between our theory and the experimental 
data for tightly entangled solutions and loosely entan- 
gled dense melts (Fig. ini- However, the Lin-Noolandi 
conjecture (and hence also our ansatz) fail for loosely en- 
tangled theta-solutions which seem to be better described 
by the theory of Colby and Rubinstein [IB, HH- An im- 
portant oversimpliflcation was recently pointed out by 
Milner [131 : by accounting only for the mean segment 
density, one implicitly treats the solvent as providing a 
"coating" which increases the effective chain diameter to 
d/Vl. As consequence, one neglects the increased prob- 
ability of entanglement formation between closely ap- 
proaching chain sections. 
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FIG. 2: Dimensionless plateau moduli G%i l^/ksT as a func- 
tion of the dimensionless number density of Kuhn segments 
PkI'k' Theory vs. experimental and scaling results as in 
Fig. [H The black solid line indicates the combination of our 
Eqs. (HI and © with = 0.06 and cq = 0.6. 



III. PRIMITIVE PATH ANALYSIS 

Are these effects preserved in an explicit primitive 
path analysis? To answer this question we have gen- 
erated and analyzed model polymer liquids correspond- 
ing to the experimental data: (i) tightly entangled so- 
lutions of zero-diameter wormHke chains (WLCs) with 
10 < pkI%: < 10^, (ii) dense melts of flexible bead-spring 
chains with $ = 0(1) and 1 < pkI]^ < 40, and (iii) 
model theta-solutions generated by ehminating a frac- 
tion 0.1 < $ < 0.925 of the chains from equihbrated (3l| 
conformations of dense bead-spring chain melts. 



A. Bead-spring polymer solutions and melts 

For the loosely entangled regime we used the model 
and procedure described in Refs. [l^. [23|. Monomers are 
modeled as spheres of diameter a interacting through 
a purely repulsive 6-12 Lennard- Jones (LJ) potential, 
which is short ranged and purely repulsive. The polymers 
are formed by connecting beads via non-linear FENE 
springs. The average bond length is 6 = 0.97cr. The 
parameter choice ensures that two chains cannot cross 
each other in dynamic simulations. Monodisperse poly- 
mer melts of Af = 80-500 chains of length 50 < < 700 
at a bead density of p = 0.85cr~'^ are studied. By in- 
troducing a small intrinsic bond bending potential, h is 
varied between 1.82a and 3.34(T, for details see Ref.[3l|. 
For dense systems, we extended the data from Ref. [23| 
by one additional data point for larger intrinsic bending 
stiffness close to the isotropic nematic transition. Fur- 
thermore, we equilibrated two Kremer-Grest theta so- 
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lutions M = 200 chains of length N = 700 at a bead 
density of p = 0.25cr~'^ and p = 0.4cr~'^. In this case, 
the LJ cutoff is set to Tc — 2.5a and simulations are car- 
ried out at ksT = S.Oeij. To obtain a large number of 
model theta-solutions whose intra- and inter-chain corre- 
lations are identical to those of the dense systems, we ran- 
domly eliminate a fraction 0.1 < <& < 0.925 of the chains 
from equilibrated [3l| conformations of dense bead-spring 
chain melts with N — 7000. The PPA is iraplemented 
into a standard Molecular Dynamics code [ll,!!^]: Chain 
ends are fixed in space, intra-chain excluded volume as 
well as bending interactions are disabled, and chain con- 
traction is induced by cooling the system toward T = 0. 



Here we have instead used Eq. ([5]) which can be written 

as 



'-'jV ''K 



2/5 
cgc^ 



8/5 



(9) 



Eqs. (HI) and ^ are equivalent, if app/lx oc ^pp/lx oc 
^/{pkIk)^ for the dense, loosely entangled melts 

where Eq. ^ was used in Ref. PPA. However, Eq. ([8]) 
clearly fails in the tightly entangled regime and we also 
found less good agreement with experimental data for es- 
timates of plateau moduli for theta-solutions (data not 
shown). 



B. Entangled solutions of zero-diameter WLCs 

Data covering the crossover to the tightly entangled 
regime were mainly obtained using Monte Carlo tech- 
niques. We have generated and analyzed semi-dilute 
solutions of infinitely thin worm-like chains with 1 < 
PkIk — Bending of the chain was penalized by the 

Hamiltonian H/ksT ~ ^(^p/b)J2i=i '^i ' Ui+i, where b 
denotes the segment length, Ip = K/ksT is the persis- 
tence length of the continuum worm-like chain with the 
bending modulus k, and Uj the unit vector along the axis 
of the «-th cylinder. Chains of this type have a Kuhn 
length of Z^V^ = 2/(1 + b/lp - cothlp/b) - 1 and can 
be efficiently generated by simple sampling. In the ab- 
sence of inter-chain correlations for zero chain diameter, 
fully equilibrated semi-dilute solutions are obtained by 
placing TV randomly generated chains into a cubic box of 
size LI^^ with periodic boundary conditions. As a general 
rule, the chains were chosen to be long enough to be mul- 
tiply entangled and the linear dimension of the simulation 
box exceeds the typical chain extensions. For the PPA 
we used a Brownian Dynamics/force-biased Monte-Carlo 
simulation code for entangled worm- like chains, which re- 
jects moves leading to chain crossing. Starting from the 
solution conformations, the positions of the chain ends 
are fixed. Chain contraction is induced by setting chain 
stiff'ness as well as the equilibrium extension of the seg- 
ments to zero and reducing the temperature. 



C. Estimating plateau moduli from the PPA 



IV. RESULTS AND DISCUSSION 



In Fig. HVl we show our PPA results for the plateau 
moduli in comparison to the experimental data. The 
PPA identifies the location of the crossovers from the 
theta-solution to the dense melts and from the loosely to 
the tightly entangled regime with no adjustable param- 
eters. In particular, we obtain quantitative agreement 
with the experimentally measured moduli by adjusting a 
single parameter (equal to cqc^^^ = 0.2) for the strength 
of the elastic response. Before we take a closer look at the 
results for loosely entangled theta solutions, we discuss 
the microscopic length scales characterizing the primitive 
path mesh. 

In Fig. Uv] we show the Kuhn lengths app and mesh 
sizes ^pp which we extracted for the various model sys- 
tems. As for the predicted moduli, the results for tightly 
entangled WLC solutions and for dense melts of fiexi- 
ble bead-spring chains are in excellent agreement with 
our theoretical results. Note that there are no adjustable 
parameters in addition to the values of — 0.06 and 
cq = 0.6 derived above from the experimentally mea- 
sured plateau moduli. 

Interestingly, the PP mesh size ^pp is essentially given 
by the packing length, p. This length scale was intro- 
duced [32| as the characteristic spatial distance where 
diflFerent polymers start to interpenetrate, i.e. the length 
scale where the intra- and interchain monomer pair cor- 
relation functions coincide: gintra{p) = ginterip)- In the 
systems we investigate interchain correlations are small. 



so that we use ginterip) = 1 or 



dL(p.lK) ^ 



Following again Refs. [1J,|23| we measure mean-square calculate the packing length. Using 
internal distances (^p (|i — j\) as a function of chemical 



distance |i — j| to determine the Kuhn length, Opp, and 
the contour length, Lpp of the primitive paths. 

For fiexible polymers the standard relation between the 
plateau modulus and the length scales characterizing the 
primitive path mesh is [l| 



L{r, Ik) 



1^(1 + L/Ik) 
ril + r/lK) 



P PkIk to 



(10) 
(11) 



the result is 



ks-i \app 



(8) 



P_ 
Ik 



i + yi + jPKt'K) 

{PKl\) 



(12) 
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FIG. 3: Dimensionless plateau moduli G% Ix/kBT as a 
function of the dimensionless number density of Kuhn seg- 
ments PkI'k- PPA vs. experimental and scaling results as 
in Fig. [H Additional symbols indicate PPA results for (i) 
tightly entangled, zero-diameter worm like chains (blue □), 
(ii) Kremer-Grest melts and theta-solutions (red and green 
filled symbols) for bead-spring polymers with intrinsic stiff- 
ness □ <A<v<O<0. Furthermore, we have included 
PPA results for model theta-solutions (open red symbols) cre- 
ated by eliminating chains from the corresponding melts. 



The limiting cases are again simple to understand. For 
Ik, the chain contour length L contained in a volume 
is much longer than the Kuhn length Ik- The chains 
follow a random walk statistics with (i?^) — IkL. In this 
case, the packing length is given by p = {pchain{R^)y^ 
[s^ l and with p/Ik oc Upp/lK oc £,pp/lK oc the 
packing length turns out to be the only relevant length 
scale in loosely entangled systems [ill, [131 ■ In the oppo- 
site limit, when p <^ Ik the single chain statistics cor- 
responds to a rigid rod with {R^) = L^. In this case, 
P = {pcha^n{RY'^)~^'^ = l/(pif Zif ) , i.e. p is identi- 
cal to the mesh size of the solution as well as the primitive 
path mesh. In contrast to the previous case, the Kuhn 
length Opp = Ik of both chains and primitive paths re- 
mains as a second, independent length scale. 

As a last point, we come back to the case of loosely 
entangled (model) theta-solutions. Both, experimental 
and PPA results for plateau moduli show clear devi- 
ations from our theory (respectively the packing con- 
jecture) G%lj^/{kBT) (X (pkIk)^ ^nd the correspond- 
ing extrapolation of the behavior of loosely entangled 
dense melts (Fig. lIVp . Qualitatively, sufficiently di- 
lute systems behave in accordance with a power law 
G% l^/{kBT) oc (pkIk)" where the exponent a is in 
the range of thepredictions of the binary contact 
model (a = 2) [TtI, [l8| and the Colby- Rubinstein model 
(a = 7/3) [2l|, [22|- However, for the relatively dense 
model theta solutions which we analyzed using the PPA, 
we see no evidence for immediate departures from the 




FIG. 4: Kuhn length app/lK and mesh size ^pp/Ik for the 
primitive path mesh as a function of the dimensionless num- 
ber density of Kuhn segments PkIr- The solid lines indicate 
the combination of Eq. ^ with Eqs. ^ and ^ respectively 
with — 0.06 and cg = 0.6. The dashed line represents 
the packing length Eq. lfT2|l . Symbols indicate PPA results 
for zero-diameter wormlike chains (blue), dense bead-spring 
melts (green) and model theta-solutions (red). 



packing line of the form G%{^) = G%{^ = 1)*", i-e. 
for a family of curves with identical slope but different 
prefactor essentially determined by the plateau moduli 
of the undiluted melt. Rather, data obtained for model 
theta solutions map fairly well onto each other and con- 
tinue to map onto the melt results as long as pkIk > ^■ 
Deviations only occur for lower concentrations. 

The comparison can be extended to the length scales 
characterizing the primitive path meshes. Inserting the 
predictions for = L^/Ik of the binary contact and 
the Colby- Rubinstein model into Eqs. ^ and ^ (which 
follow directly from the definitions and hold indepen- 
dently of the applicability of the Semenov theory to 

the primitive path mesh) yields Upp/lK oc (pkIk) 
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, ^pp/Ik oc (pkIk) and app/lx oc (prIk) 

^pp/Ik c>c {pxlji) respectively. In particular, one 

obtains for the ratio app/£,pp oc {prIkY^^ and Upp/S^pp oc 

{prI^Y^^ respectively. Thus, in the limit in question 
where (prIj^) can become arbitrarily small, these mod- 
els predict that the area mapped out by two sub- 
sequent Kuhn segments of a primitive path is on av- 
erage traversed by less than one other primitive path 
which could serve as an obstacle. The results of the ac- 
tual PPA shown in Fig. IIVI allow no definite conclusion. 
Naive extrapolation would indeed lead to Opp < ^pp in 
the limit (prI^) 0. However, we note that we have 
never observed this inversion. Interestingly, the devia- 
tions from our theory occur once p exceeds Ir and the 
results of the explicit PPA for theta-solutions are com- 
patible with a crossover to a pure packing scenario with 
p/Ir = Qpp/lR = £,pp/Ir cx 1/{prI^) where the PP di- 
rections before and after an entanglement point become 
completely uncorrelated. 

V. SUMMARY AND CONCLUSION 

There is a wide spectrum of entangled polymer liquids 
whose single chain structure is characterized by a single 
microscopic length scale, the Kuhn or persistence length. 
Our starting point was a compilation of experimentally 
measured plateau moduH for semi-dilute theta solutions 
of synthetic polymers, the corresponding dense melts 
above the glass transition or crystallization temperature, 
and solutions of semi-fiexible (bio) polymers such as 
f-actin or suspensions of rodhke viruses. Together these 
systems cover the entire range from loosely to tightly 
entangled polymers. In this article, we have demon- 
strated excellent agreement between the experimentally 
measured plateau moduli and our results derived from 
a primitive path analysis [23| of corresponding model 
polymer liquids. This is strong evidence supporting our 
working hypothesis that the PPA may be regarded as 
a tool to renormalize a loosely to a tightly entangled 
system: while the relation between the microscopic 
structure and the viscoelastic properties is different for 
the three classes of entangled polymer Hquids included in 



the present study, these differences vanish in the course 
of the PPA. In particular, experimental properties can 
be calculated by applying relations for tightly entangled 
systems to the primitive path mesh independently of the 
character of the original system (note that Eq. |[5|) for the 
plateau modulus reduces to the standard expression [l| 
Eq. ([HI) in the case of loosely entangled melts; however, 
the latter fails in the tightly entangled regime). This 
ansatz provided a new explanation of the Lin-Noolandi 
packing conjecture [l^, [l^l and allowed us to derive 
simple analytical expressions for macroscopic (Fig. U]) 
and microscopic (Fig. IIVP entanglement properties 
which are in excellent agreement with experimental 
and simulation data over a wide range of reduced 
polymer densities. Future work has to show whether the 
Colby and Rubinstein scaling [2l| in theta-solutions is a 
cross-over effect or valid asymptotically and whether it is 
possible to account for density fluctuations by combining 
the arguments and methods presented in this paper 
with the scaling analysis of loosely entangled solutions 
presented in Ref. [23|- Judging from the available 
experimental and PPA results, this might be necessary 
to "disentangle" the influence of the various crossovers 
in the experimentally relevant parameter range. 
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